PRL Preprint 



Fermionic light in common optical media 
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Recent experiments have proven that the response to short laser pulses of common optical me- 
dia, such as air or Oxygen, can be described by focusing Kerr and higher order nonlinearities of 
alternating signs. Such media support the propagation of steady solitary waves. We argue by both 
numerical and analytical computations that the low power fundamental bright solitons satisfy an 
equation of state which is similar to that of a degenerate gas of fermions at zero temperature. Con- 
sidering in particular the propagation in both O2 and air, we also find that the high power solutions 
behave like droplets of ordinary liquids. We then show how a grid of the fermionic light bubbles can 
be generated and forced to merge in a liquid droplet. This leads us to propose a set of experiments 
aimed at the production of both the fermionic and liquid phases of light, and at the demonstration 
of the transition from the former to the latter. 
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In suitable optical media, light has been argued to ac- 
quire material properties. A long dated example is the 
equivalence of the paraxial propagation of a laser pulse 
in a Kerr medium with the time evolution of a super- 
fluid Bose-Einstein condensate, due to the identity of the 
nonlinear Schrodinger equation with the Gross-Pitaevskii 
equation [1[. More recently, optical induction has been 
used to create photonic crystals [2]; a photonic system has 
been designed that may undergo a Mott insulator to su- 
perfluid quantum phase transition [3]; soliton solutions for 
light propagation in Cubic- Quintic (CQ) nonlinear media 
have been shown to behave like ordinary liquids [4|, [5j. 

On the other hand, recent experimental and theoret- 
ical works have proven that the response to ultrashort 
laser pulses of common optical media, such as air or Oxy- 
gen, can be described by focusing Kerr[6| and higher or- 
der nonlinearities of alternating signs Q, which have also 
been argued to provide the main mechanism in filament 
stabilization, instead of the plasma defocusing[7]. 

In this letter, we demonstrate by analytical and nu- 
merical computations that such media can support the 
propagation of steady solitary waves that appear in two 
clearly different phases. The low power solitons are gov- 
erned by the same equation of state than a degenerate 
gas of fermions. We will call such a system "Fermionic 
Light". On the other hand, the high power localized 
states satisfy the Young-Laplace equation that governs 
the formation of droplets in ordinary liquids, similarly to 
the result that was recently obtained for the CQ model 
We also show how to generate a grid of fermionic light 
bubbles and make it turn into a liquid light droplet. 

We will consider the (paraxial) propagation along the 
^-direction of a linearly polarized laser beam, so that the 
complex electric field component ty(x,y,z) satisfies the 
nonlinear Schrodinger equation 

' 1 ~ 9 - T - k An^ = 0, (1) 



dz 2/cono 



= d 2 /dx 2 + d 2 jdy 2 is the transverse Laplace opera- 
tor, and fco = 27r/Ao is the mean wavenumber in vacuum, 
where Ao is the central wavelength of the laser source 
which will be fixed to Ao = 800nra throughout this pa- 
per as in the experiment of Ref. [7] . For the optical media 
that have been studied in Ref. Q , the nonlinear correction 
An to the refractive index can be expanded as 



An ~ n 2 |^| 2 + n 4 |^| 4 + n 6 |^| 6 + n 8 |*| 



(2) 



with alternating sign coefficients n 2 , uq > and 724, n$ < 
0, that contribute to focusing and defocusing respec- 
tively. Taking into account that the values of the second- 
order dispersion and multiphoton-absorption coefficients 
for the air are k" = 0.2 fs 2 /cm and /? = 1.27 x 1(T 126 
cra 17 /W 9 [8], respectively, we have checked that both ef- 
fects do not lead to significant corrections in our results. 
To be concrete, we will perform most of the numerical cal- 
culations in the case of 2 as the propagation medium, 
taking the mean values obtained in the experiment , 
n 2 = 1.6 x l(T 19 cra 2 /W, n 4 = -5.2 x Kr 33 cra 4 /W 2 , 
n 6 = 4.8 x l(T 46 cra 6 /W 3 , n 8 = -2.1 x l(T 59 cra 8 /W 4 . 
For comparison, however, we will also mention the results 
that we obtain for air, using the corresponding values for 
the coefficients n 2(? that are also given in Ref.[7[. 

We will search for finite localized solutions of Eq. (pp) 
of the form &(x,y,z) = <£(r)e _w , where r 



y/x 2 



where no is the linear refractive index of the medium, 



and fx is the propagation constant. Fig. 1 shows the 
result of our numerical computation for $(0) in the exis- 
tence domain of such solitons in Oxygen (/x^ < \i < 0), 
where 11^ — —0.096 cm -1 , and for the radial profiles 
3>(r) of three of them (see left inset in FigH]), correspond- 
ing to the low power (black solid line), moderate power 
(red dashed line) and high power regimes (blue dashed- 
dotted line) respectively. We express the transverse spa- 
tial variables x, y in terms of the adimensional coordi- 
nates £, x, measured in units of (n^/no) 1 / 2 (ko^) -1 . The 
amplitude |3>(0)| is measured in units of (n 2 /n4) 1 / 2 . 
Like in the CQ case, \i can be identified with the chem- 
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FIG. 1: [Color online] Plot of the central amplitude $(0) vs. 
fi/fioo for the family of fundamental eigenstates of Eq.(pQ). 
Note that <I>(0) has a maximum for p/poo ~ 0.7, due to the 
defocusing nonlinear it ies. Left inset: radial profiles of three 
eigenfunctions of Eq.JT]) with p/poo = 0.004 (black solid), 
p/poo = 0.3 (red dashed) and p/poo = 0.97 (blue dashed- 
dotted) respectively. Right inset: radial pressure profiles cor- 
responding to these solutions. The low-power profile is mag- 
nified by a factor 10 3 for clarity. Labeled points on the main 
curve refer to the eigenstates displayed within the insets. 



ical potential of an equivalent thermodynamical two- 
dimensional system ^ of N = J pdxdy = J \ <&\ 2 dxd 
particles, described by the Landau's grand potential [1 
Q = — J pdxdy, where the pressure field p is now 



1 



2k n 



|V ± $| 2 +/i|$| : 



* $2(g+l) 

h^n 2q - ■ - (3) 
g=i 



q+1 



For our optical system, /i, TV, Q and p correspond to the 
propagation constant, the power, the Lagrangian leading 
to Eq. (pp), and the Lagrangian density, respectively. 

The right inset of FigJT] shows the pressure distribu- 
tions p(r), measured in units of kon^n^ 2 , corresponding 
to the stationary states discussed above. Notice that all 
the radial distributions display both positive and nega- 
tive pressure regions, being this a necessary condition for 
the existence of solitary waves. For the low power soliton 
(black solid line), corresponding to small values of |/i|, we 
can obtain an analytical expression by using the varia- 
tional method with the ansatz 3>(A, r) = Aexp(— r 2 /R 2 ). 
The values of A(p) and R(p) that minimize Q for a given 
value of p can then be used to compute the pressure 
distribution. Taking the first non- vanishing order in \p\ 
and inverting the dependence R(p), we find a simple an- 
alytical approximation for the central pressure p c as a 
function either of the radius R 
or of the central density p c 
the beam intensity), 

Pc = k^ 2 w = bkon2 ^ (4) 

with a = 4 and b = 1/4. This relation is similar to the 
equation of state of a degenerate gas of fermions of mass 



Y // 2(r 2 ) of the soliton, 
| $(0)| 2 (corresponding to 



m at zero temperature. In fact, if we apply the general 
definition of the Fermi momentum [To[ Pp to a two di- 
mensional system, we obtain Pp = h^J2itp 1 with p the 
density of the Fermi gas. As a consequence, the pres- 
sure, defined as the average force on a unit orthogonal 
line in the gas, can be obtained from the average kinetic 
energy as follows 



p = p(E kin ) 



P 



J Pf P 2 PdP _ vh 2 



2m 



C PdP 



2m' 



(5) 



which shows the same dependence with p 2 as Eq. (j4]). 
This proves the formal analogy of our low power solitons 
with a degenerate Fermi gas in the central region around 
r = 0, which is arbitrarily large in the limit p — >• (cor- 
responding to large radius R — » oo). For these reasons, 
we will call 'Fermionic' the phase where the pressure is 
proportional to p 2 . 
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FIG. 2: [Color online] Plot of the logarithms of p c vs. R for 
the fundamental solitons (black line). The fermionic behav- 
ior of Eq.j4j) (red dotted) and the liquid YL equation (green 
dashed) are compared with the numerics. The labeled points 
correspond to the same eigenstates displayed in Fig[T] For 
each value of R, two different outcomes are possible depend- 
ing on the power, corresponding to two different phases of 
light. Inset: plot of ln(p c ) vs. ln(p c ) (black solid) overlapped 
with Eq.JD) (red dashed) for comparison. 



Note that in the limit p 
tation gives a constant N - 



our variational compu- 
2n which is consistent 



with the known result for the power flow leading to the 
collapse threshold in a Kerr medium [9]. The magnitude 
of this power N lies in the range of few GW in both O2 
and air, and can be interpreted as the threshold for the 
existence of the Fermionic Light solitons. 

Fig |2] shows the numerical computation of p c as a func- 
tion either of R (black solid line) or p c (inset: black solid 
line) for all the nodeless solitary states of the model. The 
lower branch, corresponding to the low power solitons, is 
in excellent agreement with the dependence described by 
Eq.((4j), as it can be inferred from the fitting (red dotted) 
straight line with slope si ow = —4. Furthermore, in the 



inset of Figj2] we also show the quantitative agreement 
between theory and numerics by comparing p c with p c 
instead of R. In this case, the slope of the straight line is 
s inset — 2 thus demonstrating the quadratic dependence 
on p c given by Eq.(j4]). However, the correct numerical 
values of the constants are a — 2.5 and b = 0.29, in 
reasonable agreement with the result of the variational 
method given above. We have checked numerically that 
the asymptotic behavior represented by the red lines in 
Fig 12] is practically independent on the higher order non- 
linearities U2 q (q > 1), as predicted by the theory. In 
particular, these results can be directly applied to both 
Kerr and CQ models. 

On the other hand, the high power localized solu- 
tions exhibit top-flat profiles with an inhomogeneous neg- 
ative pressure profile on the border, similar to those 
of the solitons appearing in the CQ model[5]. As a 
consequence, the gradient term in Eq.([T]) can be ne- 
glected close to the origin, and we get p — — /coAn(O) = 
~~ Y^q=i n 2g < £(0) 2<? . By generalizing the argument of 
Ref.[5] including the higher order nonlinear it ies, we ob- 
tain that these states obey the celebrated Young-Laplace 
(YL) equation [11], p c = 2a / ' describing the behavior 
of usual liquid droplets. The value of the surface tension 



is 



n 2 



Jo 



-Poo® 2 -ko^ n 2q— | 1 ) 



$2(g+l) 



1/2 



q=l 



1 



(6) 

being p^ and the asymptotic values corresponding 
to the R — » 00 droplet, that can be computed by solving 
the equation p c = (neglecting the Laplacian term). For 
the propagation in Oxygen, we have obtained the follow- 
ing variational estimations: p^ — — 0.247(/co^2^4 X ) = 
-0.096 cm-\ $L = 0.712(n 2 n4 x ) = 2.19 x 10 13 W/cm 2 , 



a = 0.0715(n|n; 



-3/2 -1/2 
n 



) = 4.88 x 10 9 W/cm 2 . These 



results, obtained assuming a top-flat function, are in ex- 
cellent agreement with the computation given in FigJU 
In Figl21 we show that our numerical solution in the 
case of O2 satisfy the YL equation (green dashed line) 
with very good accuracy for a wide range of values of 
p. On the other hand, for the propagation in air, the 
liquid light phase would correspond to a higher intensity, 
[$^] a = 2.98 x 10 13 W/cm 2 , with [p^ = -0.116 cm" 1 
and [a] a = 7.16 x 10 9 W/cm 2 , and the YL equation would 
still be valid. 

As we have seen above, the propagation of self-guided 
light beams in media like Oxygen can occur in two clearly 
separated phases, satisfying two different equations of 
state. In fact, we have checked that qualitatively sim- 
ilar results can also be obtained for the CQ case, and 
occur whenever the nonlinear refractive index displays 
a single well-defined maximum as a function of the in- 
tensity. It would then be interesting to demonstrate the 
possibility of a transition between the fermionic bubbles 




FIG. 3: [Color online] Isosurface intensity plot of the dy- 
namical phase transition from a square grid of fermionic light 
solitons to an individual liquid light soliton in O2. All soli- 
tons in the grid are eigenstates with p/poo = 0.3 displayed in 
FigCQ They are forced to merge in the center of the compu- 
tational window by means of an external harmonic potential. 
We can see three pseudocolor plots displaying the initial con- 
dition (left, 77 = 0), the collapse of the light bubbles (middle, 
77 = 200) and the final state (right, 77 = 500). In the latter we 
observe a flat-top soliton with radius R ~ 10 and p c ~ 0.76 
arising after the massive coalescence of the fermionic soli- 
tons. The width of the square spatial domain displayed is 
w ^ x 200. 



and the liquid droplets of light. A suggestive analogy is 
that of the collapse of a star, that occurs when the grav- 
itational interaction overcomes the Fermi pressure of the 
electrons. We can obtain a qualitatively similar result in 
the case of light propagation in Oxygen, by compressing 
the fermionic bubbles using a harmonic potential, lead- 
ing to the generation of a liquid light droplet. Fig. 3 
shows the result of our simulation in O2 • The initial state 
(see left snapshot in Fig. [3]) consists of a regular grid of 
fermionic light bubbles with p/poo = 0.3 (see their ra- 
dial profile in Fig. [I]), with a separation between nearest 
neighbors A^ x = 40. We include an external harmonic 
potential V(£,x) = f (£ 2 + X 2 ) with K = 5 x 10~ 5 , in 
Eq.(pQ) in order to induce a net force acting on the grid 
with the aim of making all the fermionic solitons to col- 
lide in the center of the computational window (£, x) = 0. 

This parabolic potential can be obtained by inducing 
in the medium a " gas lens" [12| , which can be constructed 
with an electrically heated pipe through which passes a 
laminar flow of gas. By controlling the differential heat- 
ing at the boundaries and the velocity of the flow, a 
parabolic refractive index gradient can be obtained as 
documented in Ref.[12[. On the other hand, we have also 
checked numerically that the same qualitative result can 
be obtained by using a glass lens, provided that its focal 
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length / is at least ten times greater than the Rayleigh 
length of the input beams, which in the case depicted in 
Fig. 3 corresponds to / = 5 m. 
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FIG. 4: [Color online] Proposal for the generation of the initial 
condition of Fig[3] in experimental setups. The (a) pseudo- 
color plot shows the bidimensional cosine squared-like phase 
distribution to be imprinted onto the ingoing beam. We as- 
sume that the beam is much wider than the filaments, and 
restrict our simulation to a square domain £, x £ [—75, 75] 
where the wavefront is assumed to be homogeneous. In snap- 
shots (b)-(c) (77 = 100, 200) we observe how the beam profile 
starts to destabilize, leading to the formation of a grid of spa- 
tial solitons at 77 = 300 (see snapshot (d)). 



phase transition between them would be an affordable 
challenge in real experiments. 

In conclusion, we have proven that common media (O2, 
air) can support the propagation of solitary waves that 
appear in two clearly different phases with unequal phys- 
ical properties, namely the low power "Fermionic" light, 
satisfying an equation of state similar to that of a degen- 
erate gas of fermions, and the high power "Liquid" light, 
obeying the YL equation. We have then shown how a grid 
of the fermionic light bubbles can be generated and forced 
to merge in a liquid droplet. We think that the possible 
experimental validation of our proposal could also pro- 
vide an independent way to corroborate the deep change 
in the understanding of the filamentation process in gases 
that was proposed in Ref.[7]. Furthermore, these results 
in air pave the way for the improvement of recent exper- 
iments on laser-induced water condensation 14], built on 
top of these new robust light distributions. 



For simplicity, we have prepared the initial state with 
exact solutions of Eq.(pp) to reduce the excitation of lin- 
ear radiative modes. We have checked that such an ini- 
tial guess can be generated experimentally by means of 
the modulational instability of a low power probe pulse. 
In fact, FigH] shows the results of a numerical simulation 
where we have added an initial cosine squared- type phase 
distribution (see left snapshot in FigH]) to a homogeneous 
plane wave in order to excite a regular grid of solitons, as 
in Refs. [13|. As shown in snapshots b) to d), this allows 
us to control the spatial location of the optical filaments 
generated during the beam breakup due to modulational 
instability. The outgoing grid of spatial solitons is de- 
picted in snapshot d). Note its similarity with the initial 
state of FigEI 

Let us come back to the results of FigUl The massive 
coalescence of all the fermionic bubbles occurs after a 
finite propagation distance 77 = 200 (see middle snapshot 
in Fig|3]), measured in units of n^k^ri^ 2 • As a result, a 
unique filament structure with large radius arises, as it 
can be appreciated in Fig|3l We have checked that this 
soliton is a flat-top eigenstate with radial perturbations 
coming from the transition process. In fact, we have 
estimated its logarithmic radius to be around 2.3 and 
its peak density p c ~ 0.76 (measured in units of 712/^4, 
thus corresponding with an intensity 2.34 x 10 13 W/cm 2 ), 
which is clearly on the liquid light branch displayed in 
FigEI For OxygenQ, we have considered a propagation 
distance of about 13 m. Although such a distance seems 
to be huge for a real experiment, we have verified that 
by stretching the external potential this distance can be 
realistically reduced by almost one order of magnitude. 
Finally, we note that the use of pressurized O2 or air may 
also help to enhance the nonlinear optical response. For 
all these reasons, we conclude that the demonstration of 
the existence of both Fermionic and Liquid light and the 
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